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Abstract 

Using an asymptotically exact real space renormalization procedure, we find 

that the Heisenberg antiferromagnetic spin-1 chain undergoes an impurity 

driven second order phase transition from the Haldane phase to the random 

singlet phase, as the bond distribution is broadened. In the Haldane phase 

and near the critical point, there is a Griffiths region in which the gap is filled 

and the susceptibility diverges in a non-universal manner. The correlation 

length critical exponent is ~ 2.3. 
Pacs: 75.10.J, 75.30.H, 75.50.E 
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Recently, there has been tremendous interest in the antiferromagnetic (AF) spin-1 chain, 
inspired by the famous conjecture by Haldane |jl|] that integer-spin chains behave quite 
differently from half-odd-integer-spin chains. For example, in the absence of disorder, the 
spin-1 chain has short-range spin-spin correlations in the ground state and an excitation 
gap [0] whereas the spin-1/2 chain is critical. The ground state of the spin-1 chain also 
has a novel string-topological order |^. Some of these results have been experimentally 
confirmed 

Randomness is always present in real materials. Theoretical work has demonstrated that 
randomness dramatically affects the physical properties of the AF spin-1/2 chain [^,|^ and 
other random one-dimensional magnetic systems In this letter we report a systematic 

theoretical study on the effects of bond randomness on the AF spin-1 chain, based on the 
real space renormalization group scheme developed by Ma et al. (see also Ref. |]10|) and 
extended by Fisher |Q. 

Our main result is that in the presence of bond randomness, there are two distinct phases 
in the AF Heisenberg spin-1 chain, separated by a second order critical point. The nature 
of these two phases are described below. 

The ground state in the Haldane phase in the absence of randomness is well described 
by the valence bond solid (VBS) state [|ll|. In the VBS state each spin-1 is composed of 



two symetrized spin-1/2 objects. The spin-1/2 objects form singlets with spin-1/2 objects 
on neighboring sites. This state resembles the ground state of a dimerized spin-1/2 chain. 
The ground state is nondegenerate, there is an excitation gap, as well as a very stable 
topological structure. Thus we expect the Haldane phase and its topological structure to be 
stable against weak bond randomness ||] . On the other hand, when the randomness is strong 
and the distribution of bond strength is broad, the origonal spin-1 objects coupled by strong 
bonds form inert singlet pairs and generate effective further neighbor AF couplings. An 
asymptotically exact real space renormalization group (RG) analysis [^,0 shows that in this 
case the system flows toward a random singlet (RS) phase with wnwersa/ thermodynamic 
properties and power law behavior in averaged spin-spin correlations. In order to study the 



transition from the Haldane phase to the random singlet (RS) phase, we have extended 
this RG scheme so that it may be used in both phases. We find the transition between 
these two phases is second order. The extended RG scheme becomes asymptotically exact 
in the low energy limit at the critical point, as well as in the RS phase. Thus we are 
able to extract exact information about the critical point. For examples, we find as the 
randomness strength approaches the critical point from the Haldane phase, the average 
spin-spin correlation length diverges in a power law manner with exponent z/ = 2.3. 
The string-topological order parameter vanishes with a power law exponent 7 = 2i/ 4.6. 
Consider the Hamiltonian 

H = Y^JSi-^i+i, (1) 

i 

where Si are spin-S operators and Jj are random coupling constants (assumed positive unless 
otherwise specified). Let us begin with the case that randomness is strong and the width 
of the distribution of J is very broad (on a logarithmic scale). In this case we can use 
the decimation renormalization group procedure developed by Ma, Dasgupta, and Hu |^ 
(MDH) for the special case of S* = 1/2. We first pick the largest bond in the system, say 
J2 between spins 2 and 3. Since this is such a strong bond, and the width of the bond 
distribution is broad, bonds Ji and J3 will likely be much weaker than J2- Thus to lowest 
order in Ji and J3, spins 2 and 3 form a singlet pair and become unimportant at low energies 
(on scales much smaller than J2). The major physical effect of the existence of spins 2 and 
3 is to generate an induced coupling between neighboring spins 1 and 4. Hi_^ = J14S1 • S4 
where J14 = |S'(S' + l)JiJzl J2 + 0(1/J|). This formula is correct even if Ji and/or J3 are 
ferromagnetic (F) as long as their magnitude is much less than J2. The generated bond 
is typically much weaker than all three original bonds. Thus the effect of this decimation 
procedure is to get rid of the strongest bond (and also its two neighbors) in the system, 
generate a weaker bond between the spins neighboring the decimated ones, and lower the 
overall energy scale. Fisher has shown rigorously that as one proceeds with the RG 
scheme the width of bond distribution broadens and the accuracy of the procedure improves 
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as the energy scale is lowered. It becomes asymptotically exact in the long-distance, low- 
energy limit where the bond distribution flows toward an energy scale dependent stable fixed 
point distribution. Fisher names the phase characterized by this universal bond distribution 
the random singlet (RS) phase. Thus in the presence of strong enough bond randomness, 
a spin-S" Heisenberg AF chain will be in the RS phase, irrespective of the size of S and 
independent of S being an integer or a half-odd-integer. 

The above real space RG procedure may not be accurate when the initial bond distribu- 
tion is not very broad. In this case there is a finite probability that Ji and/or J3 is of similar 
strength as J2, in which case the lowest order perturbation is not sufficient. What one needs 
to do in this case is to pick the segment in the chain in which all spins are coupled by strong 
bonds (while the segment itself is only weakly coupled to the rest of the chain), solve the 
spectrum of the segment, and keep only the low energy states, and consider their coupling 
to the rest of the chain. In the case of a spin-1/2 chain, the ground state for a segment is 
either a singlet (for even segments) or a doublet (for odd segments), which are separated 
from higher energy states by a gap. In the former case the segment will be inert, and merely 
mediates a weak effective AF coupling between the two spins neighboring the segment, while 
in the latter case the segment may be modeled by an effective spin-1/2 at low energy which is 
coupled to the rest of the chain antiferromagnetically. Thus the structure of the RG scheme 
remains the same as the MDH procedure, even if greater care is necessary in the beginning. 
It is therefore believed that the RS phase correctly describes the long-distance, low-energy 
physics of the spin-1/2 chain, even if there is only weak randomness 0. 

The situation is very different in the case of the spin-1 chain. For a finite segment of 
the spin-1 chain with no disorder, there are two effective half spins localized near the two 
edges of the segment, and the coupling between them is g{l) ~ ( — l)'Ja', where / is the 
length of the segment and a < 1 |]12[. Thus the coupling between the two half spins in the 



same segment may be either ferromagnetic or antiferromagnetic, and decays exponentially 
with the length of the segment. The coupling of an edge spin to the rest of the chain, 
however, remains antiferromagnetic. The low energy states of a segment are the singlet 



and triplet states formed by these two edge spins (and both of them must be kept); other 
states will be separated from them by the Haldane gap of order J . We expect this scenario 
to remain correct even when there is weak randomness in the bond strength in such finite 
segments. Thus even though the spin chain is composed of spin-1 objects, the effective 
degrees of freedom at low energies are actually half spins. This is also true when there is 



no randomness. The valence bond solid state is also composed of half spins [Tl||. We may 
describe the low energy physics of a random spin-1 chain using an an effective Hamiltonian 
in terms of these half spins, with the following structure: the chain consists of half spins only] 
the even bonds are taken from an antiferromagnetic bond distribution, and odd bonds are 
taken from a distribution containing antiferromagnetic and ferromagnetic bonds. Physically 
the even bonds are couplings between edge half spins of neighboring segments (which are 
always AF), and odd bonds are the coupling between edge spin in the same segment. Such 
a effective description is particularly accurate for the special case of a spin-1 chain with a 
bimodal bond distribution: most of the bonds are of strength J, while a small but finite 
fraction of bonds have a much smaller strength J'. In this case the system can indeed be 
viewed as a collection of weakly coupled segments of uniform chains. The couplings between 
edge half spins are random because the length of the segments are random. We believe this 
model correctly describes the long-distance, low-energy physics of a random spin-1 chain. 
In particular, the original spin-1 Hamiltonian may be recovered by setting all odd bonds 
strongly ferromagnetic |l^ . 



To study this model, we develop an extended version of the MDH real space renormal- 
ization scheme, which properly accounts for strong ferromagnetic bonds. In our scheme, at 
any stage of RG, the energy scale VL is set by the strongest antiferromagnetic bond in the 
system. All even bonds, which are all AF, will be weaker than Vl. We separate the odd 
bonds, which can be either F or AF, into two groups: group A consists of all AF bonds and 
those F bonds that are weaker than fi, while group B consists of F bonds that are stronger 
than VL. The new spin decimation procedure works in the following way. Find the strongest 
AF bond in the system, say Jj. Hi is odd, then its neighbors are both weaker AF bonds. 



We just follow the usual MDH decimation procedure. If i is even and both neighbors belong 
to group A we can again follow the MDH decimation procedure. If one of the bonds, say 
Jj+i, belongs to group B, we solve the 3-spin cluster problem of Si, Si+i and Si+2, and keep 
the low energy states, which are a doublet. The doublet may be modeled by a new spin-1/2, 
which couples to the rest of the chain. If both neighbors of Jj belong to group B, we solve 
the 4-spin cluster including spin i — 1, i, i + 1 and i + 2. The ground state is a singlet, with 
an excitation gap of order Jj. Thus it drops out at low energy and mediates an effective 
AF coupling between spins i — 2 and i + 3. It is easy to show that this procedure keeps the 
original structure of the system; i.e., even bonds AF and odd bonds F or AF. 

The flow equations for the distributions of bonds in this approximate RG scheme are 

dPe{J,n) 



+ [Pe{n, Q)(i - A^'(Q)) - Po{n, n)] p,{j) 
-6{n- j)Pe{fi,n), 



(2) 



d\l J-u 



- 6{n - j)Poin, n) + 2Pe{n, n)N{n)Po{-j, n) 



+ 



Po{n,n)-p,{n,a){i-N\a)) p,{j), 



(3) 



dN{n) 

dn 



Po{n,n)-Pe{^,^){'i--N'^{^)) N{n) 



+ Po{-n,n). 



(4) 



Here Pe{J,fl) is the (normalized) probability distribution of even bonds with fl > J > 0, 
Po{J, ^) is the probability distribution of odd bonds with Q, > J > — Q, and N{Q) is the 
fraction of odd bonds that are strongly ferromagnetic J < —Vt. Po{J,Q) and N[Q) are 
related by the normalization condition J^QdJPo{J,fl) + N{fl) — 1. Anticipating that the 
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bond distribution will become very broad in the low energy limit, we have neglected factors 
of order 1 in the strength of generated bonds. They become irrelevant in the asymptotic 
limit p. We have also assumed that the ferromagnetic bonds that are stronger than Q are 
much stronger than Q, so that two spin-1/2 objects connected by a strong ferromagnetic 
bond form a spin-1 object. Again this assumption is valid in the asymptotic limit (fi ^ 0), 
and simplifies the solution of clusters including strong F bonds. 

The density, p{fl), of spins that have not yet paired into singlets at scale fl satisfies 

dp{n) 



[Pein, n) [i + N^n)] + Po{n, n)]p{n). (5) 



These spins are essentially free at temperatures that are higher than Q. All thermodynamic 
quantities can be determined from p{^l) |]5|, p!4| , p!5 . 



It is natural to combine the distributions in the following way 

p_(j, n) = PoiJ, n) - Po{-j, ^]), (6) 

so that Q+{J, fi) is normalized. We find that fixed point solutions to the flow equations are 
power laws Pe{x) = Pe~^^ and Q+{x) = Qe'^^ where x = log{Q/J) and P_ = |n,T3 . 
The variables P, Q, N, and p obey 

dT 2 ^ 

^ = -il-N)QP 

dT ^ '^2 ' 

^ = -[(i + iv^)p + i^g]p, (7) 

where F = log(fij/fi) (fij is the initial cutoff of AF bonds). We find there are two classes 
of stable fixed points, corresponding to two stable phases. They are the random singlet 
phase (P = V^^,Q = Qo,N = l,p oc F~^), and the Haldane phase (P = Po,Q = 0,N = 
0,p QC Q^°, where < Pq < 1 is a nonuniversal number). There is also an unstable fixed 
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point (near which there is a single relevant operator): {P = Q = 2T~^,N = |,p oc T~^). It 
describes the critical point. In the following we describes the physical nature of these fixed 
points. 

In the random singlet phase, all the odd bonds become F bonds much stronger than the 
AF even bonds. All the spin-l/2's are ferromagnetically combined into spin-l's. The spin- 
I's then couple into singlets over all length scales. The ground state and thermodynamic 
properties are the same as for the spin-1/2 random singlet state studied previously The 
disorder averaged spin-spin correlation function C{r) decays as r~^, and the susceptibility 
take the universal form in the low temperature limit: x ^ [Tlog^T]"^. There is no gap 
because the bond distribution has weight at J = 0. 

In the Haldane phase all odd bonds (F and AF) become much weaker than the even 
bonds, only spin-l/2's remain in the system and they form singlets only over even bonds. 
The system may be viewed as a set of uncoupled dimers. The spin-spin correlations decay 
exponentially with a finite correlation length. There is also long range string-topological 
order p[. This phase is analogous to the random dimer phase in random dimerized AF 
spin-1/2 chains ^ and the ground state resembles the valence bond solid state. The flow 
equations describe the Griffiths region of the Haldane phase where there is no gap and the 
susceptibility diverges as a power law with a nonuniversal exponent x ~ T~^^~^°\ The flow 
equations are only valid when the disorder is broad, and so can not describe the crossover 
from gapped to gapless behavior within the Haldane phase as the randomness is increased. 

To determine critical exponents we consider small perturbations near the unstable fixed 
point: 

p = l(i + 6,r') 

Ar = l(l + 5„r^) (8) 

and expand the fiow equations to linear order in 6. There are two irrelevant perturbations 
(A = — 1,A = ^^-j^) and one relevant perturbation (A+ = relevant fiows, if 
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5n > 0, the odd bonds are stronger than the even bonds, the density of spin-l's is increasing, 

and the system flows to the random singlet fixed point. If 5„ < the even bonds are stronger 

than the odd bonds, the density of spin-l's is decreasing, and the system flows to the random 

Haldane phase. The crossover from critical to Haldane behavior occurs at the energy scale 

i_ 

where ST^ 1. The energy scale at which this occurs is Fq = 5 ^+ . The density of spins 
at this scale is po = Tq so the correlation length is ^ ~ pg ^ = 6^'^ where u = ^ 2.3. 
Following the same analysis as in Ref. P], the string-topological order parameter defined 
as i 



T = lim ( ^0 



exp 








i<k<j 





^0 , (9) 



scales as T oc {—6n)^'^, in the Haldane phase near the critical point. Just like the RS 
fixed point for the spin-1/2 chain 0, the bond distributions become infinitely broad on a 
logarithmic scale. Thus our scheme becomes asymptotically exact at the critical point, and 
the critical exponents we find here are exact. 

Before closing we briefly discuss relations between the present work and some existing 
papers. Our model is related, but different from that studied by Westerberg et al. |]^. In 
their model, ferromagnetic bonds may appear in both even and odd bonds, and spins of 
arbitrarily large size appear at low energy. The special topological structure of our model 



prevents this from happening in the present system. Boechat et al. [0 also anticipated the 
existence of a random singlet phase for strong randomness. As in our previous work 0, 
they also find spontaneously dimerized chains are unstable against weak randomness. They 
did not, however, address the weak randomness regime of the spin-1 chain, and the phase 
transition that we discuss here. The correlation length exponent v ^ 2.3 we find here is 
extremely close to that of the delocalization transition in integer quantum Hall systems. Lee 
ir?! showed that the integer quantum Hall transition may be mapped onto the dimerization 
transition in the pure SU(0) spin chain. It is unclear at this stage, whether this is merely 
a coincidence, or there exists a fundamental physical reason, that these two apparently 
different transitions have the same critical exponent. 



To summarize, in this letter we determined the critical properties of the randomness 
driven phase transition of the spin-1 chain. For weak randomness, the spin-1 chain is in the 
Haldane phase and the ground state resembles the valence bond solid state. The ground 
state has topological order, spin-spin correlations decays exponentially and there is a gap to 
the excited states. For broader distributions the gap is filled in and the correlation length 
increases but the topological order persists. Eventually the critical point is reached and the 
correlation length diverges. The correlation length exponent is u — ~ 2.3. Beyond 

the critical point the topological order is zero, the disorder averaged spin-spin correlations 
decay algebraically, and the ground state resembles the random singlet state. 
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